**************************************
**************************************
* Replication file for Lacina, B and Lee, C 2013
* �Culture Clash or Democratic Peace?: 
* Results of a Survey Experiment on the Effect of Religious Culture and Regime Type on Foreign Policy Opinion Formation,�
* Foreign Policy Analysis
**************************************
**************************************

*LOAD FEELING THERMOMETER DATA

use "LacinaLee_FPA_replicationFig1and2.dta", clear

notes

*GENERATE FIGURE 1

twoway (lfit thermometer polity) (scatter thermometer polity, sort msymbol(none) mlabp(0) ///
 mlabel(abbrev)), ytitle(Average feeling thermometer rating (1-100)) ylabel(,nogrid) ///
 xtitle(`label') legend(off) scheme(lean2) 
 
*GENERATE CORRELATIONS REPORTED IN FIGURE 1
 
reg thermometer polity
 
pcorr thermometer polity christian muslim

*GENERATE FIGURE 2

twoway (lfit thermometer christian) (scatter thermometer christian, sort msymbol(none) mlabp(0) ///
 mlabel(abbrev)), ytitle(Average feeling thermometer rating (1-100)) ylabel(,nogrid) ///
 xtitle(`label') legend(off) scheme(lean2) 

twoway (lfit thermometer muslim) (scatter thermometer muslim, sort msymbol(none) mlabp(0) ///
 mlabel(abbrev)), ytitle(Average feeling thermometer rating (1-100)) ylabel(,nogrid) ///
 xtitle(`label') legend(off) scheme(lean2) 
 
*GENERATE CORRELATIONS REPORTED IN FIGURE 2

reg thermometer christian
reg thermometer muslim

pcorr thermometer polity christian muslim

*LOAD EXPERIMENTAL RESULTS

use "LacinaLee_FPA_replicationdata.dta", clear

notes

*Table 1

tab threatp

tab trustp

tab threatt

tab trustt

*Data shown in Figure 3

tab threatp if demop==1

tab threatp if demop==2

tab threatp if religp==1

tab threatp if religp==2

*Data shown in Figure 4

tab threatt if demot==1

tab threatt if demot==2

tab threatt if religt==1

tab threatt if religt==2

*Table 2

tab threatp if demop==1

tab threatp if demop==2

tab threatp if religp==1 & demop==1

tab threatp if religp==1 & demop==2

tab threatp if religp==1

tab threatp if religp==2 & demop==1

tab threatp if religp==2 & demop==2

tab threatp if religp==2

ktau threatp religp

ranksum threatp, by(religp) porder

ktau threatp religp if demop==1

ktau threatp religp if demop==2

ktau threatp demop

ranksum threatp, by(demop) porder

ktau threatp demop if religp==1

ktau threatp demop if religp==2

gen notchristiandemp=1

replace notchristiandemp=0 if demop==1 & religp==1

ktau threatp notchristiandemp

ranksum threatp, by(notchristiandemp) porder

*Table 3

tab threatt if demot==1

tab threatt if demot==2

tab threatt if religt==1 & demot==1

tab threatt if religt==1 & demot==2

tab threatt if religt==1

tab threatt if religt==2 & demot==1

tab threatt if religt==2 & demot==2

tab threatt if religt==2

ktau threatt religt

ranksum threatt, by(religt) porder

ktau threatt religt if demot==1

ktau threatt religt if demot==2

ktau threatt demot

ranksum threatt, by(demot) porder

ktau threatt demot if religt==1

ktau threatt demot if religt==2

gen notchristiandemt=1

replace notchristiandemt=0 if demot==1 & religt==1

ktau threatt notchristiandemt

ranksum threatt, by(notchristiandemt) porder

*Data shown in Figure 5

tab trustp if demop==1

tab trustp if demop==2

tab trustp if religp==1

tab trustp if religp==2

*Data shown in Figure 6

tab trustt if demot==1

tab trustt if demot==2

tab trustt if religt==1

tab trustt if religt==2

*Table 4

tab trustp if demop==1

tab trustp if demop==2

tab trustp if religp==1 & demop==1

tab trustp if religp==1 & demop==2

tab trustp if religp==1

tab trustp if religp==2 & demop==1

tab trustp if religp==2 & demop==2

tab trustp if religp==2

ktau trustp religp

ranksum trustp, by(religp) porder

ktau trustp religp if demop==1

ktau trustp religp if demop==2

ktau trustp demop

ranksum trustp, by(demop) porder

ktau trustp demop if religp==1

ktau trustp demop if religp==2

gen notchristiandemp=1

replace notchristiandemp=0 if demop==1 & religp==1

ktau trustp notchristiandemp

ranksum trustp, by(notchristiandemp) porder

*Table 5

tab trustt if demot==1

tab trustt if demot==2

tab trustt if religt==1 & demot==1

tab trustt if religt==1 & demot==2

tab trustt if religt==1

tab trustt if religt==2 & demot==1

tab trustt if religt==2 & demot==2

tab trustt if religt==2

ktau trustt religt

ranksum trustt, by(religt) porder

ktau trustt religt if demot==1

ktau trustt religt if demot==2

ktau trustt demot

ranksum trustt, by(demot) porder

ktau trustt demot if religt==1

ktau trustt demot if religt==2

gen notchristiandemt=1

replace notchristiandemt=0 if demot==1 & religt==1

ktau trustt notchristiandemt

ranksum trustt, by(notchristiandemt) porder


